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ABSTRACT 

We present results for neutron star models constructed with a new equation of state for 
nuclear matter at zero temperature. The ground state is computed using the Auxiliary 
Field Diffusion Monte Carlo (AFDMC) technique, with nucleons interacting via a semi- 
phenomenological Hamiltonian including a realistic two-body interaction. The effect 
of many-body forces is included by means of additional density-dependent terms in 
the Hamiltonian. In this letter we compare the properties of the resulting neutron-star 
models with those obtained using other nuclear Hamiltonians, focusing on the relations 
between mass and radius, and between the gravitational mass and the baryon number. 
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1 INTRODUCTION 

While real neutron stars are very complicated objects, their 
main global properties can usually be well-approximated 
by considering simple idealized models consisting of perfect 
fluid in hydrostatic equilibrium. If rotation can be neglected 
to a first approximation (as is the case for the spin rates 
of most currently-known pulsars) then the model can be 
taken to be spherical and its structure obtained by solv- 
ing the Tolman-Oppenheimer-Volkoff (TOV) equations, en- 
abling one to calculate, for example, the stellar mass as a 
function of radius or of central density. However, this re- 
quires specification of an equation of state (EOS) for the 
neutron-star matter which clearly plays a fundamental role 
in determining the properties of the resulting models. The 
EOS needs to take account of the strong spin-isospin cor- 
relations induced by reali stic interacti ons, with particular 
regard to the tensor ones (|Raffeltll 19961 ). 

The EOS can in principle be computed by means of 
many-body theories using effective density-dependent in- 
teractions, such as those given by Skyrme forces but the 
phenomenological nature of these can be a disadvantage 
in making reliable calculations of neutron-star properties 
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jRikovska Stone et al.l [20031 ). Making a microscopic calcu- 
lation based on a Hamiltonian describing the properties of 
light nuclei and symmetric nuclear matter (SNM), is both 
challenging and of great relevance. Variational approaches 
using correlated basis functions (C BF) and Fermi Hy- 
per N etted Chain techniques (FHNC) (|Fantoni fc Fabrocinil 
Il998h are good candidates for doing this but the strong 
spin-isospin dependence of the nuclear Hamiltonian requires 
the introduction of approximations into the FHNC scheme 
which cannot be fully controlled, such as the si ngle operator 
chain (SOC) (|Pandharipande fc Wiringalll979l ). 

The A uxiliary Field Diffusion M onte Carlo method 
(AFDMC) (|Schmidt fc Fantonil [l999l ) has been employed 
for computing the ground state energy of nuclei, giving 
very good agreement with ot her existing accurate tech- 
niques in few-body physics (|Gandolfi et al.l l2007al ) and 
highlighting important limitations of other many-body 
theories used for cal culations of both nu clear structure 
( Gandolfi eral]|2007bl ) and neutron matter (|Gandolfi et al.l 
l2009l ). In this letter we present results from an AFDMC 
calculation of the EOS of /3-equilibrium matter (relevant 
for neutron stars), based on a new class of non-relativistic 
Hamiltonians using nuclear potentials with a realistic two- 
body interaction, taken from the Urbana-Argonne scheme 
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l|Wiringa et alJ Il995l ) and incorporating many-body inter- 
actions via density-dependent terms. 

Modern nuclear Hamiltonians are usually constrained 
so as to reproduce proper ties of light nucle i as measured 
in laboratory experiments (|Pieper et al.ll200ll 'l. but the den- 
sities involved there are much lower than those found in 
neutron-star cores for which it becomes necessary to in- 
troduce a three-bo dy potential, such as UIX or IL1-IL5 
jPieper et al.l 1200 ih . However, these give very large and 
very different energy contributions to the EOS of pure 
neutron matter (PNM) at high de nsity (jSarsa et al.ll2003l : 
iFantoni et alll2007l: iGandolfil |2007| '1. Also, recent AFDMC 
results bv lGandolfi et all (|2007bf )indicate that SNM is not 
well reproduced by the FHNC/SOC techniques used to con- 
strain these Hamiltonians, which are therefore not very suit- 
able for calculating neutron-star models. 

Because of this, we have proceeded here in a different 
way, using density-dependent terms to simulate many-body 
forces with the forms used coming from explicit integration 
of these forces over the variables of particle 3 for the three- 
body force, over those of particles 3 and 4 for the four-body 
force, and so on. In bulk matter, the neutron and proton 
densities p„ and Pp can be taken to be constant quantities, 
whereas in confined systems, like nuclei, they are operators. 

As a first step in this direction, we have re- 
visited the ol d three-parameter den sit y-dep endent 
LP model of Lagaris fc Pandharipandg (|l98ll ) and 



[Friedman fc Pandharipandd ( 198ll ), with the values of 
the free parameters fixed so as to reproduce the saturation 
point and compressibility of SNM. The LP calculations 
were performed using the FHNC/SOC approximation 
and we then re-fitted the density-dependent term using 
the AFDMC calculations. By considering the chemical 
potentials, we have constructed an EOS for a mixture of 
protons, electrons and muons in /3-equilibrium and have 
then used the resulting EOS for constructing neutron-star 
models. We stress that (i) the ground states of nuclear 
matter as a function of pp and pn have been calculated with 
the same AFDMC method which has been shown to be 
very accurate for calculating the ground states of various 
nucleonic systems, and (ii) the many-body part of the 
nuclear interaction has been self-consistently determined by 
these solutions. 



roperties and to give t he best fit to NN scattering data 

Wiringa fc Piepej|20o3 ) . 

The many-body interactions are represented by density- 
dependent factors of the structural form given by the LP 
model. The resulting potential, denoted as DD6', is given 
by the following six two-body components 



TNA(p) 



372 p e 




(1) 



with 71 , 72 and 73 being fixed so as to reproduce the exper- 
imental values of the saturation density po = 0.16 fm~^, the 
binding energy per particle Eq = —16 MeV and the com- 
pressibility K = 9pg {d'^E{p)/dp'^) ^ 240 MeV. (Note 
that in this paper we follow the nuclear- physics convention 
of using p to refer to a particle number density rather than 
a mass density.) 

Many nuclear matter calculations of the 1980s were 
made with the LP model interaction, and gave very good 
agreement with data for large nuc lei, as discussed by 
iBenhar fc Pandharipandd (|l993l ) and iPandharipande et al.l 
( 1997^ . Afterwards, mainlv because of the difficulty of treat- 
ing the density as an operator when dealing with nuclei, this 
model was forgotten and attention was switched to the de- 
velopment of increasingly sophisticated three-body poten- 
tials. On the basis of the results obtained with AFDMC 
quantum simulations, we believe that we should instead now 
proceed in the direction of constructing increasingly sophis- 
ticated density-dependent two-body effective interactions, 
which may also account for N-body interactions with A'^ ^ 3. 

A particular limitation of DD6' is the exclusion of the 
non-local components of AV18. The most important of these 
is the spin-orbit term, particularly when dealing with nuclei. 
We estimate that the missing non-local components con- 
tribute to the EOS o f neutron matter by no more than 5% 
(|Gandolfi et al.ll2009l ). Work towards introducing the spin- 
orbit components is in progress. A second important point 
concerns the structural form of the density-dependent terms 
and one may want to increase the number of Feynman dia- 
grams included in the construction of UIX and contributing 
to the many-body interaction. 



2 THE MODEL 

We model the EOS by simulating nuclear matter using a 
finite number of interacting nucleons in a periodic box. The 
number of particles is chosen from among the magic num- 
bers giving a rotationally invariant wave function for the 
corresponding non-interacting system and the volume of the 
box is fixed by the density. For the two-body interaction, w e 
take the Argonne AV6' potential (|Wiringa fc Piepej[200l ). 
which includes the four central spin-isospin components and 
the two tensor ones. The six components of the long range 
OPEP potential are fully included, whereas only the first six 
of the 18 components of the i ntermediate and sho rt range 
parts of the AV18 interaction (|Wiringa et al.|[l995l ). Vj{rij) 
and Vg{rij), are kept. The corresponding amplitudes Ip and 
Sp are re-fitted so as to correctly reproduce the deuteron 



3 CALCULATION OF THE EQUATION OF 
STATE OF NUCLEAR MATTER 

We computed the ground state of the system using the 
AFDMC method (for details, see lGandolfill2007l and refer- 
ences therein), simulating SN M with A = 28 nucleon s in a 
periodic box, as described bv iGandolfi et all (|2007bt ). The 
next magic number of nucleo ns providing small fin ite-size 
corrections is 132, as shown bv lGandolfi et al] |20o3), which 
requires an unjustified computational effort given the re- 
striction to local components in the DD6' model. Table [1] 
gives the values of the free parameters of DD6' correspond- 
ing to the best fit to experimental data for po, Eo and K, and 
compares these with the values for the original LP model. 

The AFDMC density-dependent terms give more at- 
traction than in the original LP model, consistent with the 
fact that FHNC/SOC overbinds SNM. The TNA term is 
~ 30% larger at po and more than the double at 5po, giving 
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Table 1. The AFDMC results for the free parameters of the DD6' 
interaction as compared with the original values of the LP model 
calculated within the FHNC/SOC approximation. 



parameter 


FHNC/SOC 


AFDMC 


71 


0.15 fm^ 


0.10 fm^ 


72 


-700 fra^ 


-750 fm<5 


73 


13.6 fm^ 


13.9 fm^ 



more attraction, and exp(7ip) — 1 is '--^ 30% smaller over the 
whole range (po,5po), giving less repulsion. 

We find that the AFDMC results for the binding energy 
per nucleon of SNM at densities larger than ~ 0.08 fm~^ can 
be very well described by: 

Esnm{p) = Eo + a{p - pof + h{p - po)'e^(''-'"'' , (2) 

where Eo = -16.0 MeV, po = 0.16 fm-^ a = 520.0 MeV 
fm^ 6 = -1297.4 MeV fm^ and 7 = -2.213 fm^ This 
parametrization was chosen to represent the EOS of nuclear 
matter, reproducing properties constrained b y terrestrial ex- 
periments on nuclei (|Danielewicz et al]|2002l '). 

The DD6' Hamiltonian was then used to compute the 
EOS of PNM, by making a simulation with 66 neutrons in 
a periodic box. The EOS for nuclear matter as a function of 
the proton fraction Xp = pp/ p is then parametrized as 

E{p, x-p) = Esnm{p) + Cs (^j^ (1 - 2xj,f . (3) 

The two extra parameters of the symmetry energy term, Cs 
and 7s, were obtained by fitting E{p, Xp — 0) to the AFDMC 
resuh for PNM. This gives Cs = 31.3 MeV and 7^ = 0.64. 
Typical values for these parameters have be en quoted as 
Cs ~ 31 - 33 MeV and 7^ « 0.55 - 0.69 by IShettv et al.1 
ll2007l^ and a s Cs = 31.6 MeV and 7^ « 0.69 - 1.05 by 
IWorlev et al.1 (|2008l '). It should be noted that usually the 
symmetry energy is constrained over a range of densities 
typical of nuclei, whereas we have here fitted the parame- 
ters over a very wide density range. This means that the 
parametrization of eq. ^ should be accurate up to very 
high densities. 

In high-density matter, neutrons can produce protons 
and electrons by /3 decay and so the equilibrium configura- 
tion can have a non-zero proton/neutron ratio, modifying 
the EOS away from that for PNM. The equilibrium concen- 
tration of protons Xp can be computed by imposing 

= Mp + Me I (4) 

where fii is the chemical potential (for neutrons, protons 
and electrons respectively). For doing this, we consider the 
electrons as comprising a relativistic Fermi-gas: 

fi. = [mi + h\3np,f^']'^\ (5) 

and charge neutrality imposes that pe = Xpp, where p is 
the total nucleon density. The chemical potentials of the 
neutrons and protons are derived from equation (|3]), and 
equation Q is then solved to give Xp as a function of p. 
Another consideration is that when /ie becomes larger than 
the muon mass, the production of muons is favoured, and 
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Figure 1. The equation of state for symmetric nuclear mat- 
ter (dashed line), pure neutron matter (dot-dashed line) and /3- 
equilibrium nuclear matter with both electrons and muons (full 
line) and with electrons only (dotted line). The points show the 
AFDMC results for symmetric nuclear matter (SNM, squares) 
and pure neutron matter (PNM, circles) which have been used to 
fit the EOS. In the inset, we show the proton fraction Xp plotted 
as a function of the total nucleon density, computed consider- 
ing the presence of just electrons as negatively-charged particles 
(dotted line) and with both electrons and muons (full line). 



their contribution must then also be considered. For our 
EOS, muons begin to appear at p ~ 0.133 fm~^. 

The EOSs calculated for SNM, PNM and /3-equilibrium 
matter are shown in Fig. [U with the energy per nucleon 
being plotted against the nucleon number density p; in the 
inset, we show the proton fraction Xp plotted as a function 
of p, computed considering both electrons and muons (full 
line) and only electrons (dotted hue) . The EOS for PNM is 
softer than tha t calcu lated with the AVS'-fUIX potential by 
iGandolfi etaP (j2009l ) using the same AFDMC many-body 
method and this then produces a similar behaviour for the 
/3-equilibrium matter. The main reason for this comes from 
the different treatment of many-body effects. 

At very high densities, the chemical potential of the 
nucleonic matter becomes larger than the threshold for cre- 
ation of heavier particles. Such states are due to the up and 
down quarks transforming to strange quarks, so that parti- 
cles with strangeness (hyperons) start to appear. A realistic 
EOS should include these when they appear and this can 
seriously modify the structure of the star. We do not in- 
clude this in our present calculations (although we make 
some comment about its likely effects in the next section); 
we leave inclusion of this until a subsequent paper. 



4 RESULTING NEUTRON STAR MODELS 

When the EOS of the neutron-star matter has been speci- 
fied, the structure of an idealized spherically-symmetric neu- 
tron star model can be calculated by integrating the Tolman- 
Oppenheimer-Volkoff (TOV) equations: 
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dP 

dr 



G[m{r) + 4Tvr^P/c^][e + P/c^] 



r[r ~ 2Gm(r)/c2] 



dr 



(6) 
(7) 



where P = p {dE /dp) and e 



p{E + mjv) are the pressure 
and the energy density, mjv is the average nucleon mass, 
rn(r) is the gravitational mass enclosed within a radius r, 
and G is the gravitational constant. The solution of the TOV 
equations for a given central density gives the profiles of p, 
e and P as functions of radius r, and also the total radius 
R and mass M = m{R). A sequence of models can be gen- 
erated by specifying a succession of values for the central 
density. In Fig. [2] we plot the mass M (measured in so- 
lar masses Mq) as a function of the radius R (measured in 
km), as obtained from calculations with four different pre- 
scriptions for the EOS: the /3-equilibrium and PNM EOSs 
discussed in Sections 2 and 3, the equivalent one for PNM 
with just two-body interactions (using AV6'), and a previous 
one from ICandolfi et all (|2009l '). for PNM with three-body 
interactions (using AV8'+UIX). Models to the right of the 
maximum of each curve are stable to radial perturbations 
and these are the ones of interest for astrophysical neutron 
stars. The maximum mass obtained with the two-body in- 
teraction AV8' is very similar to that for AV6'. 

It is interesting to make a comparison between these 
curves so as to see the changes caused by introduction of the 
various different features. The solid curve (/^-equilibrium) is 
our best proposal for the neutron-star EOS but it can be seen 
that it differs only very little from the pure neutron matter 
EOS (where the radii for a given mass are just slightly larger 
within the main range of interest). There is a considerable 
difference, however, with respect to the previous AV8'-fUIX 
curve for pure neutron matter, with the maximum mass be- 
ing reduced from ~ 2.5 M@ to ~ 2.2 Mq and the radii in 
the main region of interest also being substantially reduced. 
This reflects the effective softening of the EOS caused by 
the different treatment of many-body effects. 

An objective of this type of work is to attempt to 
constrain microphysical models for neutron-star matter by 
making comparison with astronomical observations (see 
iLattimer fc PrakashI [20071 ') . This is just starting to be pos- 
sible now and further progress is anticipated within the 
next few years. At present, the only neutron stars for which 
masses are accurately known are the components of the best- 
observed binary pulsars, for which timing measurements give 
results correct to many significant figures. The maximum 
mass for any of these is the 1.44 1 Mq for the Hulse/Taylo r 
binary pulsar PSR 19134-16 (see IWeisberg fc Tavloj|2005h . 
However, there is accumulating evidence for higher masses, 
particularly for neutr on stars in binary systems together 
with white dwarfs (see lRansom et al]|2005l ) and there is now 
a widespread belief that the maximum should probably be 
in the range 1.8 Mq — 2.1 Mq (at least when the rotation 
is sufficiently slow, as is the case for almost all pulsars so 
far observed). At high enough densities, it is expected that 
the composition of the matter would change because of the 
appearance of either hyperons or deconfined quarks, both 
of which are likely to decreas e the maximum mass (see, for 
example, Sc 
lAkmal et al. 



ulze et al ] |2006l for the case of hyperons, and 
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Figure 2. Predicted neutron-star masses (in units of Mq) plot- 
ted as a function of stellar radius (in km). Four different equa- 
tions of state are considered: those discussed in this paper for 
/3-equilibrium matter (full line) and pure neutron matter (PNM, 
DD6', dashed line), the equivalent one for PNM with just two- 
body interactio ns (PNM, AV6' , dot-d ash-dashed line), and a pre- 
vious one from [Gandolfi et al ] liooi), for PNM with three-body 
interactions (PNM, AV8'-|-UIX, dot-dashed line). 



mass for /3-equilibrium matter is pc ~ 1.2 fm~^ which is well 
within the range where these changes are likely to have oc- 
curred, and so we expect that the maximum mass would be 
slightly lower than that shown in Fig. [51 This brings us well 
within the expected range. 

Observational constraints for the radius are more prob- 
lematic, but one of the best of these seems to be the in di- 
rect constraint suggested by [Podsiadlowski et al.l (|2005l ) in 
the case of the less massive component of the double pulsar 
PSR J0737-3059. If, as seems likely, this neutron star was 
the product of an electron-capture supernova, then the to- 
tal pre-coUapse baryon number of the stellar core is rather 
precisely known from model calculations and, since only a 
very small loss of material is expected to have occurred in 
the subsequent collapse, the baryon number of the neutron 
star is itself also well-known. Together with the very accu- 
rate value for the gravitational mass (calculated from pulsar 
timing) , this can be used to place a quite stringent constraint 
on the EOS. The baryon number A of a neutron-star model 
can be readily calculated from 



(8) 



matter). The central density corresponding to our maximum 



dAjr) 2 / 2Gm\-2 

— =47rpr 1 —] 

dr \ rc-' I 

which needs to be solved together with equations ((Gjl and 
((7| (we recall that p is here the baryon number density). 
In practice, it is convenient to talk in terms of the fearyonic 
mass, defined as Mo = mjvj4(J?), rather than A(IV) itself: 
the difference between the baryonic mass and the gravita- 
tional mass depends on the compactness of the neutron star, 
and hence indirectly on the radius. If M is plotted against 
Mo for a given EOS, then the curve needs to pass through 
a certain error box in order to be consistent with the ob- 
servations for PSR J0737-3059, subject to the assumptions 
being made in the analysis. This plot is shown in Fig. [3] for 
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Figure 3. The gravitational mass of stellar models plotted as a 
function of the baryonic mass. See the caption of Fig.[2]for details. 
The shaded rectangle is the error box inferred from observations 
of the lower-mass component of the doubl e pulsar PSR J0737- 
3059 and associated modelling (following iPodsiadlowski et al. 
I2OO5I. with the modification mentioned in the text) . iKitaureTetTL 
||2009|) have proposed a smaller error box just to the left of this, 
but we prefer the one shown here. 



the same EOSs shown in Fig. (2] (We use here the error box 
as suggested by Podsiadlowski et al. 2005, extended to the 
left by 3 X 10^^ M0 to include a plausible upper limit for 
the matter lost during the core collapse.) The curve for our 
/3-equilibrium EOS just touches the top left-hand corner of 
the error box. However, the central density of our model cor- 
responding to the mass of the pulsar concerned {M = 1.249 
M0) is more than three times nuclear matter density, by 
which point hyperons have probably already appeared, giv- 
ing some softening of the EOS which would move the curve 
slightly downwards. 



5 CONCLUSIONS 

In this letter we have presented a new equation of state for 
neutron-star matter based on microscopic calculations made 
with the Auxiliary Field Diffusion Monte Carlo technique, 
using a semi-phenomenological Hamiltonian including a re- 
alistic two-body interaction and with the effect of many- 
body forces being included by means of additional density- 
dependent terms. We have presented results from stellar 
model calculations using the new EOS and some related 
variants, focusing on the mass/radius relation and the rela- 
tion between the gravitational mass and the baryon number, 
and we have compared them with observational constraints. 



ACKNOWLEDGMENTS 

We thank S. Reddy, R. De Pietri and J.R. Stone for use- 
ful discussions. This work was supported by CompStar, a 
Research Networking Programme of the European Science 



Foundation; KES thanks SISSA for kind hospitality and ac- 
knowledges support from NSF grant PHY-0757703. Calcu- 
lations were performed using the HPC facility "BEN" at 
ECT* in Trento, under a grant for Supercomputing Projects, 
and using the HPC facility of SISSA/Democritos in Trieste. 



REFERENCES 

Akmal A., Pandharipande V. R., Ravenhall D. G., 1998, 

Phys. Rev. C, 58, 1804 
Benhar O., Pandharipande V. R., 1993, Rev. Mod. Phys., 

65, 817 

Danielewicz P., Lacey R., Lynch W. G, 2002, Science, 298, 
1592 

Fantoni S., Fabrocini A., 1998, Lecture Notes in Physics, 
510, 119 

Fantoni S., Gandolfi S., Pederiva F., Schmidt K. E., 2007, 
Recent Progress in Many-Body Theories, 11, 23 

Friedman B., Pandharipande V., 1981, Nucl. Phys. A, 361, 
502 

Gandolfi S., 2007, arXiv:0712.1364 [nucl-th] 

Gandolfi S., lUarionov A. Y., Schmidt K. E., Pederiva F., 

Fantoni S., 2009, Phys. Rev. C, 79, 054005 
Gandolfi S., Pederiva F., Fantoni S., Schmidt K. E., 2007a, 

Phys. Rev. Lett., 99, 022507 
Gandolfi S., Pederiva F., Fantoni S., Schmidt K. E., 2007b, 

Phys. Rev. Lett., 98, 102503 
Kitaura F. S., Janka H.-T., Hillebrandt W., 2009, Astron. 

Astrophys., 450, 345 
Lagaris L E., Pandharipande V. R., 1981, Nucl. Phys. A, 

359, 331 

Lattimer J. M., Prakash M., 2007, Physics Reports, 442, 
109 

Pandharipande V. R., Sick L, deWitt Huberts P. K. A., 

1997, Rev. Mod. Phys., 69, 961 
Pandharipande V. R., Wiringa R. B., 1979, Rev. Mod. 

Phys., 51, 821 

Pieper S. C, Pandharipande V. R., Wiringa R. B., Carlson 
J., 2001, Phys. Rev. C, 64, 014001 

Podsiadlowski P., Dewi J. D. M., Lesaffre P., Miller J. C, 
Newton W. G, Stone J. R., 2005, Mon. Not. R. Astron. 
Soc, 361, 1243 

Raffelt G. G, 1996, The stars as laboratories of fundamen- 
tal physics. University of Chicago, Chicago & London 

Ransom S. M., Hessels J. W. T., Stairs L H., Freire P. C. C, 
Camilo F., Kaspi V. M., Kaplan D. L., 2005, Science, 307, 
892 

Rikovska Stone J., Miller J. C, Koncewicz R., Stevenson 
P. D., Strayer M. R., 2003, Phys. Rev. C, 68, 034324 

Sarsa A., Fantoni S., Schmidt K. E., Pederiva F., 2003, 
Phys. Rev. C, 68, 024308 

Schmidt K. E., Fantoni S., 1999, Phys. Lett. B, 446, 99 

Schulze H.-J., Polls A., Ramos A., Vidafia L, 2006, Phys. 
Rev. C, 73 

Shetty D. V., Yennello S. J., Souliotis G. A., 2007, Phys. 

Rev. C, 76, 024606 
Weisberg J. M., Taylor J. H., 2005, in Rasio F. A., Stairs 

L H., eds, Binary radio pulsars, ASP Conference Series, 

Vol. 328 p. 25 

Wiringa R. B., Pieper S. C, 2002, Phys. Rev. Lett., 89, 
182501 



6 S. Gandolfi, A. Yu. Illarionov, S. Fantoni, J.C. Miller, F. Pederiva, K.E. Schmidt 

Wiringa R. B., Stoks V. G. J., SchiaviUa R., 1995, Phys. 
Rev. C, 51, 38 

Worley A., Krastev R G., Li B.-A., 2008, Astrophys. J., 
685, 390 



This paper has been typeset from a T^^X/ WT^^ file prepared 
by the author. 



